Comparative analysis of Wolbachia maternal transmission and localization in host ovaries.

Many insects and other animals carry microbial endosymbionts that influence their reproduction and fitness. These relationships only persist if endosymbionts are reliably transmitted from one host generation to the next. Wolbachia are maternally transmitted endosymbionts found in most insect species, but transmission rates can vary across environments. Maternal transmission of wMel Wolbachia depends on temperature in natural Drosophila melanogaster hosts and in transinfected Aedes aegypti, where wMel is used to block pathogens that cause human disease. In D. melanogaster, wMel transmission declines in the cold as Wolbachia become less abundant in host ovaries and at the posterior pole plasm (the site of germline formation) in mature oocytes. Here, we assess how temperature affects maternal transmission and underlying patterns of Wolbachia localization across 10 Wolbachia strains diverged up to 50 million years-including strains closely related to wMel-and their natural Drosophila hosts. Many Wolbachia maintain high transmission rates across temperatures, despite highly variable (and sometimes low) levels of Wolbachia in the ovaries and at the developing germline in late-stage oocytes. Identifying strains like closely related wMel-like Wolbachia with stable transmission across variable environmental conditions may improve the efficacy of Wolbachia-based biocontrol efforts as they expand into globally diverse environments.

Aedes albopictus is present in the lowlands of southern Zambia.

Identifying the current geographic range of disease vectors is a critical first step towards determining effective mechanisms for controlling and potentially eradicating them. This is particularly true given that historical vector ranges may expand due to changing climates and human activity. The Aedes subgenus Stegomyia contains over 100 species, and among them, Ae. aegypti and Ae. albopictus mosquitoes represent the largest concern for public health, spreading dengue, chikungunya, and zika viruses. While Ae. aegypti has been observed in the country of Zambia for decades, Ae. albopictus has not. In 2015 we sampled four urban and three rural areas in Zambia for Aedes species. Using DNA barcoding, we confirmed the presence of immature and adult Ae. albopictus at two sites: Siavonga and Livingstone. These genotypes seem most closely related to specimens previously collected in Mozambique based on mtDNA barcoding. We resampled Siavonga and Livingstone sites in 2019, again observing immature and adult Ae. albopictus at both sites. Relative Ae. albopictus frequencies were similar between sites, with the exception of immature life stages, which were higher in Siavonga than in Livingstone in 2019. While Ae. albopictus frequencies did not vary through time in Livingstone, both immature and adult frequencies increased through time in Siavonga. This report serves to document the presence of Ae. albopictus in Zambia, which will contribute to understanding the potential public health implications of this disease vector in southern Africa.

Rapid turnover of pathogen-blocking Wolbachia and their incompatibility loci.

At least half of all insect species carry maternally inherited Wolbachia alphaproteobacteria, making Wolbachia the most common endosymbionts in nature. Wolbachia spread to high frequencies is often due to cytoplasmic incompatibility (CI), a Wolbachia-induced sperm modification that kills embryos without Wolbachia. Several CI-causing Wolbachia variants, including wMel from Drosophila melanogaster, also block viruses. Establishing pathogen-blocking wMel in natural Aedes aegypti mosquito populations has reduced dengue disease incidence, with one study reporting about 85% reduction when wMel frequency is high. However, wMel transinfection establishment is challenging in many environments, highlighting the importance of identifying CI-causing Wolbachia variants that stably persist in diverse hosts and habitats. We demonstrate that wMel-like variants have naturally established in widely distributed holometabolous dipteran and hymenopteran insects that diverged approximately 350 million years ago, with wMel variants spreading rapidly among these hosts over only the last 100,000 years. Wolbachia genomes contain prophages that encode CI-causing operons (cifs). These cifs move among Wolbachia genomes - with and without prophages - even more rapidly than Wolbachia move among insect hosts. Our results shed light on how rapid host switching and horizontal gene transfer contribute to Wolbachia and cif diversity in nature. The diverse wMel variants we report here from hosts present in different climates offer many new options for broadening Wolbachia-based biocontrol of diseases and pests.

A Simulation Framework for Modeling the Within-Patient Evolutionary Dynamics of SARS-CoV-2.

The global impact of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has led to considerable interest in detecting novel beneficial mutations and other genomic changes that may signal the development of variants of concern (VOCs). The ability to accurately detect these changes within individual patient samples is important in enabling early detection of VOCs. Such genomic scans for rarely acting positive selection are best performed via comparison of empirical data with simulated data wherein commonly acting evolutionary factors, including mutation and recombination, reproductive and infection dynamics, and purifying and background selection, can be carefully accounted for and parameterized. Although there has been work to quantify these factors in SARS-CoV-2, they have yet to be integrated into a baseline model describing intrahost evolutionary dynamics. To construct such a baseline model, we develop a simulation framework that enables one to establish expectations for underlying levels and patterns of patient-level variation. By varying eight key parameters, we evaluated 12,096 different model-parameter combinations and compared them with existing empirical data. Of these, 592 models (∼5%) were plausible based on the resulting mean expected number of segregating variants. These plausible models shared several commonalities shedding light on intrahost SARS-CoV-2 evolutionary dynamics: severe infection bottlenecks, low levels of reproductive skew, and a distribution of fitness effects skewed toward strongly deleterious mutations. We also describe important areas of model uncertainty and highlight additional sequence data that may help to further refine a baseline model. This study lays the groundwork for the improved analysis of existing and future SARS-CoV-2 within-patient data.

Aedes albopictus is present in the lowlands of southern Zambia.

Identifying the current geographic range of disease vectors is a critical first step towards determining effective mechanisms for controlling and potentially eradicating them. This is particularly true given that historical vector ranges may expand due to changing climates and human activity. The Aedes subgenus Stegomyia contains over 100 species, and among them, Ae. aegypti and Ae. albopictus mosquitoes represent the largest concern for public health, spreading dengue, chikungunya, and Zika viruses. While Ae. aegypti has been observed in the country of Zambia for decades, Ae. albopictus has not. In 2015 we sampled four urban and two rural areas in Zambia for Aedes species. Using DNA barcoding, we confirmed the presence of immature and adult Ae. albopictus at two rural sites: Siavonga and Livingstone. These genotypes seem most closely related to specimens previously collected in Mozambique based on CO1 sequence from mtDNA. We resampled Siavonga and Livingstone sites in 2019, again observing immature and adult Ae. albopictus at both sites. Relative Ae. albopictus frequencies were similar between sites, with the exception of immature life stages, which were higher in Siavonga than in Livingstone in 2019. While Ae. albopictus frequencies did not vary through time in Livingstone, both immature and adult frequencies increased through time in Siavonga. This report serves to document the presence of Ae. albopictus in Zambia, which will contribute to the process of determining the potential public health implications of this disease vector in Central Africa.

A simulation framework for modeling the within-patient evolutionary dynamics of SARS-CoV-2.

The global impact of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has led to considerable interest in detecting novel beneficial mutations and other genomic changes that may signal the development of variants of concern (VOCs). The ability to accurately detect these changes within individual patient samples is important in enabling early detection of VOCs. Such genomic scans for positive selection are best performed via comparison of empirical data to simulated data wherein evolutionary factors, including mutation and recombination rates, reproductive and infection dynamics, and purifying and background selection, can be carefully accounted for and parameterized. While there has been work to quantify these factors in SARS-CoV-2, they have yet to be integrated into a baseline model describing intra-host evolutionary dynamics. To construct such a baseline model, we develop a simulation framework that enables one to establish expectations for underlying levels and patterns of patient-level variation. By varying eight key parameters, we evaluated 12,096 different model-parameter combinations and compared them to existing empirical data. Of these, 592 models (~5%) were plausible based on the resulting mean expected number of segregating variants. These plausible models shared several commonalities shedding light on intra-host SARS-CoV-2 evolutionary dynamics: severe infection bottlenecks, low levels of reproductive skew, and a distribution of fitness effects skewed towards strongly deleterious mutations. We also describe important areas of model uncertainty and highlight additional sequence data that may help to further refine a baseline model. This study lays the groundwork for the improved analysis of existing and future SARS-CoV-2 within-patient data.

Developing an appropriate evolutionary baseline model for the study of SARS-CoV-2 patient samples.

Over the past 3 years, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has spread through human populations in several waves, resulting in a global health crisis. In response, genomic surveillance efforts have proliferated in the hopes of tracking and anticipating the evolution of this virus, resulting in millions of patient isolates now being available in public databases. Yet, while there is a tremendous focus on identifying newly emerging adaptive viral variants, this quantification is far from trivial. Specifically, multiple co-occurring and interacting evolutionary processes are constantly in operation and must be jointly considered and modeled in order to perform accurate inference. We here outline critical individual components of such an evolutionary baseline model-mutation rates, recombination rates, the distribution of fitness effects, infection dynamics, and compartmentalization-and describe the current state of knowledge pertaining to the related parameters of each in SARS-CoV-2. We close with a series of recommendations for future clinical sampling, model construction, and statistical analysis.

Distinct Wolbachia localization patterns in oocytes of diverse host species reveal multiple strategies of maternal transmission.

A broad array of endosymbionts radiate through host populations via vertical transmission, yet much remains unknown concerning the cellular basis, diversity, and routes underlying this transmission strategy. Here, we address these issues, by examining the cellular distributions of Wolbachia strains that diverged up to 50 million years ago in the oocytes of 18 divergent Drosophila species. This analysis revealed 3 Wolbachia distribution patterns: (1) a tight clustering at the posterior pole plasm (the site of germline formation); (2) a concentration at the posterior pole plasm, but with a significant bacteria population distributed throughout the oocyte; and (3) a distribution throughout the oocyte, with none or very few located at the posterior pole plasm. Examination of this latter class indicates Wolbachia accesses the posterior pole plasm during the interval between late oogenesis and the blastoderm formation. We also find that 1 Wolbachia strain in this class concentrates in the posterior somatic follicle cells that encompass the pole plasm of the developing oocyte. In contrast, strains in which Wolbachia concentrate at the posterior pole plasm generally exhibit no or few Wolbachia in the follicle cells associated with the pole plasm. Taken together, these studies suggest that for some Drosophila species, Wolbachia invade the germline from neighboring somatic follicle cells. Phylogenomic analysis indicates that closely related Wolbachia strains tend to exhibit similar patterns of posterior localization, suggesting that specific localization strategies are a function of Wolbachia-associated factors. Previous studies revealed that endosymbionts rely on 1 of 2 distinct routes of vertical transmission: continuous maintenance in the germline (germline-to-germline) or a more circuitous route via the soma (germline-to-soma-to-germline). Here, we provide compelling evidence that Wolbachia strains infecting Drosophila species maintain the diverse arrays of cellular mechanisms necessary for both of these distinct transmission routes. This characteristic may account for its ability to infect and spread globally through a vast range of host insect species.

A male-killing Wolbachia endosymbiont is concealed by another endosymbiont and a nuclear suppressor.

Bacteria that live inside the cells of insect hosts (endosymbionts) can alter the reproduction of their hosts, including the killing of male offspring (male killing, MK). MK has only been described in a few insects, but this may reflect challenges in detecting MK rather than its rarity. Here, we identify MK Wolbachia at a low frequency (around 4%) in natural populations of Drosophila pseudotakahashii. MK Wolbachia had a stable density and maternal transmission during laboratory culture, but the MK phenotype which manifested mainly at the larval stage was lost rapidly. MK Wolbachia occurred alongside a second Wolbachia strain expressing a different reproductive manipulation, cytoplasmic incompatibility (CI). A genomic analysis highlighted Wolbachia regions diverged between the 2 strains involving 17 genes, and homologs of the wmk and cif genes implicated in MK and CI were identified in the Wolbachia assembly. Doubly infected males induced CI with uninfected females but not females singly infected with CI-causing Wolbachia. A rapidly spreading dominant nuclear suppressor genetic element affecting MK was identified through backcrossing and subsequent analysis with ddRAD SNPs of the D. pseudotakahashii genome. These findings highlight the complexity of nuclear and microbial components affecting MK endosymbiont detection and dynamics in populations and the challenges of making connections between endosymbionts and the host phenotypes affected by them.

Patterns of Population Structure and Introgression Among Recently Differentiated Drosophila melanogaster Populations.

Despite a century of genetic analysis, the evolutionary processes that have generated the patterns of exceptional genetic and phenotypic variation in the model organism Drosophila melanogaster remains poorly understood. In particular, how genetic variation is partitioned within its putative ancestral range in Southern Africa remains unresolved. Here, we study patterns of population genetic structure, admixture, and the spatial structuring of candidate incompatibility alleles across a global sample, including 223 new accessions, predominantly from remote regions in Southern Africa. We identify nine major ancestries, six that primarily occur in Africa and one that has not been previously described. We find evidence for both contemporary and historical admixture between ancestries, with admixture rates varying both within and between continents. For example, while previous work has highlighted an admixture zone between broadly defined African and European ancestries in the Caribbean and southeastern USA, we identify West African ancestry as the most likely African contributor. Moreover, loci showing the strongest signal of introgression between West Africa and the Caribbean/southeastern USA include several genes relating to neurological development and male courtship behavior, in line with previous work showing shared mating behaviors between these regions. Finally, while we hypothesized that potential incompatibility loci may contribute to population genetic structure across the range of D. melanogaster; these loci are, on average, not highly differentiated between ancestries. This work contributes to our understanding of the evolutionary history of a key model system, and provides insight into the partitioning of diversity across its range.

cifB-transcript levels largely explain cytoplasmic incompatibility variation across divergent Wolbachia.

Divergent hosts often associate with intracellular microbes that influence their fitness. Maternally transmitted Wolbachia bacteria are the most common of these endosymbionts, due largely to cytoplasmic incompatibility (CI) that kills uninfected embryos fertilized by Wolbachia-infected males. Closely related infections in females rescue CI, providing a relative fitness advantage that drives Wolbachia to high frequencies. One prophage-associated gene (cifA) governs rescue, and two contribute to CI (cifA and cifB), but CI strength ranges from very strong to very weak for unknown reasons. Here, we investigate CI-strength variation and its mechanistic underpinnings in a phylogenetic context across 20 million years (MY) of Wolbachia evolution in Drosophila hosts diverged up to 50 MY. These Wolbachia encode diverse Cif proteins (100% to 7.4% pairwise similarity), and AlphaFold structural analyses suggest that CifB sequence similarities do not predict structural similarities. We demonstrate that cifB-transcript levels in testes explain CI strength across all but two focal systems. Despite phylogenetic discordance among cifs and the bulk of the Wolbachia genome, closely related Wolbachia tend to cause similar CI strengths and transcribe cifB at similar levels. This indicates that other non-cif regions of the Wolbachia genome modulate cif-transcript levels. CI strength also increases with the length of the host's larval life stage, presumably due to prolonged cif action. Our findings reveal that cifB-transcript levels largely explain CI strength, while highlighting other covariates. Elucidating CI's mechanism contributes to our understanding of Wolbachia spread in natural systems and to improving the efficacy of CI-based biocontrol of arboviruses and agricultural pests globally.

Temperature effects on cellular host-microbe interactions explain continent-wide endosymbiont prevalence.

Endosymbioses influence host physiology, reproduction, and fitness, but these relationships require efficient microbe transmission between host generations to persist. Maternally transmitted Wolbachia are the most common known endosymbionts,1 but their frequencies vary widely within and among host populations for unknown reasons.2,3 Here, we integrate genomic, cellular, and phenotypic analyses with mathematical models to provide an unexpectedly simple explanation for global wMel Wolbachia prevalence in Drosophila melanogaster. Cooling temperatures decrease wMel cellular abundance at a key stage of host oogenesis, producing temperature-dependent variation in maternal transmission that plausibly explains latitudinal clines of wMel frequencies on multiple continents. wMel sampled from a temperate climate targets the germline more efficiently in the cold than a recently differentiated tropical variant (∼2,200 years ago), indicative of rapid wMel adaptation to climate. Genomic analyses identify a very narrow list of wMel alleles-most notably, a derived stop codon in the major Wolbachia surface protein WspB-that underlie thermal sensitivity of cellular Wolbachia abundance and covary with temperature globally. Decoupling temperate wMel and host genomes further reduces transmission in the cold, a pattern that is characteristic of host-microbe co-adaptation to a temperate climate. Complex interactions among Wolbachia, hosts, and the environment (GxGxE) mediate wMel cellular abundance and maternal transmission, implicating temperature as a key determinant of Wolbachia spread and equilibrium frequencies, in conjunction with Wolbachia effects on host fitness and reproduction.4,5 Our results motivate the strategic use of locally selected wMel variants for Wolbachia-based biocontrol efforts, which protect millions of individuals from arboviruses that cause human disease.6.

Male Age and Wolbachia Dynamics: Investigating How Fast and Why Bacterial Densities and Cytoplasmic Incompatibility Strengths Vary.

Endosymbionts can influence host reproduction and fitness to favor their maternal transmission. For example, endosymbiotic Wolbachia bacteria often cause cytoplasmic incompatibility (CI) that kills uninfected embryos fertilized by Wolbachia-modified sperm. Infected females can rescue CI, providing them a relative fitness advantage. Wolbachia-induced CI strength varies widely and tends to decrease as host males age. Since strong CI drives Wolbachia to high equilibrium frequencies, understanding how fast and why CI strength declines with male age is crucial to explaining age-dependent CI's influence on Wolbachia prevalence. Here, we investigate if Wolbachia densities and/or CI gene (cif) expression covary with CI-strength variation and explore covariates of age-dependent Wolbachia-density variation in two classic CI systems. wRi CI strength decreases slowly with Drosophila simulans male age (6%/day), but wMel CI strength decreases very rapidly (19%/day), yielding statistically insignificant CI after only 3 days of Drosophila melanogaster adult emergence. Wolbachia densities and cif expression in testes decrease as wRi-infected males age, but both surprisingly increase as wMel-infected males age, and CI strength declines. We then tested if phage lysis, Octomom copy number (which impacts wMel density), or host immune expression covary with age-dependent wMel densities. Only host immune expression correlated with density. Together, our results identify how fast CI strength declines with male age in two model systems and reveal unique relationships between male age, Wolbachia densities, cif expression, and host immunity. We discuss new hypotheses about the basis of age-dependent CI strength and its contributions to Wolbachia prevalence. IMPORTANCEWolbachia bacteria are the most common animal-associated endosymbionts due in large part to their manipulation of host reproduction. Many Wolbachia cause cytoplasmic incompatibility (CI) that kills uninfected host eggs. Infected eggs are protected from CI, favoring Wolbachia spread in natural systems and in transinfected mosquito populations where vector-control groups use strong CI to maintain pathogen-blocking Wolbachia at high frequencies for biocontrol of arboviruses. CI strength varies considerably in nature and declines as males age for unknown reasons. Here, we determine that CI strength weakens at different rates with age in two model symbioses. Wolbachia density and CI gene expression covary with wRi-induced CI strength in Drosophila simulans, but neither explain rapidly declining wMel-induced CI in aging D. melanogaster males. Patterns of host immune gene expression suggest a candidate mechanism behind age-dependent wMel densities. These findings inform how age-dependent CI may contribute to Wolbachia prevalence in natural systems and potentially in transinfected systems.

A single mutation weakens symbiont-induced reproductive manipulation through reductions in deubiquitylation efficiency.

Animals interact with microbes that affect their performance and fitness, including endosymbionts that reside inside their cells. Maternally transmitted Wolbachia bacteria are the most common known endosymbionts, in large part because of their manipulation of host reproduction. For example, many Wolbachia cause cytoplasmic incompatibility (CI) that reduces host embryonic viability when Wolbachia-modified sperm fertilize uninfected eggs. Operons termed cifs control CI, and a single factor (cifA) rescues it, providing Wolbachia-infected females a fitness advantage. Despite CI's prevalence in nature, theory indicates that natural selection does not act to maintain CI, which varies widely in strength. Here, we investigate the genetic and functional basis of CI-strength variation observed among sister Wolbachia that infect Drosophila melanogaster subgroup hosts. We cloned, Sanger sequenced, and expressed cif repertoires from weak CI-causing wYak in Drosophila yakuba, revealing mutations suspected to weaken CI relative to model wMel in D. melanogaster A single valine-to-leucine mutation within the deubiquitylating (DUB) domain of the wYak cifB homolog (cidB) ablates a CI-like phenotype in yeast. The same mutation reduces both DUB efficiency in vitro and transgenic CI strength in the fly, each by about twofold. Our results map hypomorphic transgenic CI to reduced DUB activity and indicate that deubiquitylation is central to CI induction in cid systems. We also characterize effects of other genetic variation distinguishing wMel-like cifs Importantly, CI strength determines Wolbachia prevalence in natural systems and directly influences the efficacy of Wolbachia biocontrol strategies in transinfected mosquito systems. These approaches rely on strong CI to reduce human disease.

Wolbachia in the spittlebug Prosapia ignipectus: Variable infection frequencies, but no apparent effect on host reproductive isolation.

Animals serve as hosts for complex communities of microorganisms, including endosymbionts that live inside their cells. Wolbachia bacteria are perhaps the most common endosymbionts, manipulating host reproduction to propagate. Many Wolbachia cause cytoplasmic incompatibility (CI), which results in reduced egg hatch when uninfected females mate with infected males. Wolbachia that cause intense CI spread to high and relatively stable frequencies, while strains that cause weak or no CI tend to persist at intermediate, often variable, frequencies. Wolbachia could also contribute to host reproductive isolation (RI), although current support for such contributions is limited to a few systems. To test for Wolbachia frequency variation and effects on host RI, we sampled several local Prosapia ignipectus (Fitch) (Hemiptera: Cercopidae) spittlebug populations in the northeastern United States over two years, including closely juxtaposed Maine populations with different monomorphic color forms, "black" and "lined." We discovered a group-B Wolbachia (wPig) infecting P. ignipectus that diverged from group-A Wolbachia-like model wMel and wRi strains in Drosophila-6 to 46 MYA. Populations of the sister species Prosapia bicincta (Say) from Hawaii and Florida are uninfected, suggesting that P. ignipectus acquired wPig after their initial divergence. wPig frequencies were generally high and variable among sites and between years. While phenotyping wPig effects on host reproduction is not currently feasible, the wPig genome contains three divergent sets of CI loci, consistent with high wPig frequencies. Finally, Maine monomorphic black and monomorphic lined populations of P. ignipectus share both wPig and mtDNA haplotypes, implying no apparent effect of wPig on the maintenance of this morphological contact zone. We hypothesize P. ignipectus acquired wPig horizontally as observed for many Drosophila species, and that significant CI and variable transmission produce high but variable wPig frequencies.

Pervasive effects of Wolbachia on host activity.

Heritable symbionts have diverse effects on the physiology, reproduction and fitness of their hosts. Maternally transmitted Wolbachia are one of the most common endosymbionts in nature, infecting about half of all insect species. We test the hypothesis that Wolbachia alter host behaviour by assessing the effects of 14 different Wolbachia strains on the locomotor activity of nine Drosophila host species. We find that Wolbachia alter the activity of six different host genotypes, including all hosts in our assay infected with wRi-like Wolbachia strains (wRi, wSuz and wAur), which have rapidly spread among Drosophila species in about the last 14 000 years. While Wolbachia effects on host activity were common, the direction of these effects varied unpredictably and sometimes depended on host sex. We hypothesize that the prominent effects of wRi-like Wolbachia may be explained by patterns of Wolbachia titre and localization within host somatic tissues, particularly in the central nervous system. Our findings support the view that Wolbachia have wide-ranging effects on host behaviour. The fitness consequences of these behavioural modifications are important for understanding the evolution of host-symbiont interactions, including how Wolbachia spread within host populations.

A phylogeny for the Drosophila montium species group: A model clade for comparative analyses.

The Drosophila montium species group is a clade of 94 named species, closely related to the model species D. melanogaster. The montium species group is distributed over a broad geographic range throughout Asia, Africa, and Australasia. Species of this group possess a wide range of morphologies, mating behaviors, and endosymbiont associations, making this clade useful for comparative analyses. We use genomic data from 42 available species to estimate the phylogeny and relative divergence times within the montium species group, and its relative divergence time from D. melanogaster. To assess the robustness of our phylogenetic inferences, we use 3 non-overlapping sets of 20 single-copy coding sequences and analyze all 60 genes with both Bayesian and maximum likelihood methods. Our analyses support monophyly of the group. Apart from the uncertain placement of a single species, D. baimaii, our analyses also support the monophyly of all seven subgroups proposed within the montium group. Our phylograms and relative chronograms provide a highly resolved species tree, with discordance restricted to estimates of relatively short branches deep in the tree. In contrast, age estimates for the montium crown group, relative to its divergence from D. melanogaster, depend critically on prior assumptions concerning variation in rates of molecular evolution across branches, and hence have not been reliably determined. We discuss methodological issues that limit phylogenetic resolution - even when complete genome sequences are available - as well as the utility of the current phylogeny for understanding the evolutionary and biogeographic history of this clade.

Comparative studies on speciation: 30 years since Coyne and Orr.

Understanding the processes of population divergence and speciation remains a core question in evolutionary biology. For nearly a hundred years evolutionary geneticists have characterized reproductive isolation (RI) mechanisms and specific barriers to gene flow required for species formation. The seminal work of Coyne and Orr provided the first comprehensive comparative analysis of speciation. By combining phylogenetic hypotheses and species range data with estimates of genetic divergence and multiple mechanisms of RI across Drosophila, Coyne and Orr's influential meta-analyses answered fundamental questions and motivated new analyses that continue to push the field forward today. Now 30 years later, we revisit the five questions addressed by Coyne and Orr, identifying results that remain well supported and others that seem less robust with new data. We then consider the future of speciation research, with emphasis on areas where novel methods and data motivate potential progress. While the literature remains biased towards Drosophila and other model systems, we are enthusiastic about the future of the field.

Pervasive Effects of Wolbachia on Host Temperature Preference.

Heritable symbionts can modify a range of ecologically important host traits, including behavior. About half of all insect species are infected with maternally transmitted Wolbachia, a bacterial endosymbiont known to alter host reproduction, nutrient acquisition, and virus susceptibility. Here, we broadly test the hypothesis that Wolbachia modifies host behavior by assessing the effects of eight different Wolbachia strains on the temperature preference of six Drosophila melanogaster subgroup species. Four of the seven host genotypes infected with A-group Wolbachia strains (wRi in Drosophila simulans, wHa in D. simulans, wSh in Drosophila sechellia, and wTei in Drosophila teissieri) prefer significantly cooler temperatures relative to uninfected genotypes. Contrastingly, when infected with divergent B-group wMau, Drosophila mauritiana prefers a warmer temperature. For most strains, changes to host temperature preference do not alter Wolbachia titer. However, males infected with wSh and wTei tend to experience an increase in titer when shifted to a cooler temperature for 24 h, suggesting that Wolbachia-induced changes to host behavior may promote bacterial replication. Our results indicate that Wolbachia modifications to host temperature preference are likely widespread, which has important implications for insect thermoregulation and physiology. Understanding the fitness consequences of these Wolbachia effects is crucial for predicting evolutionary outcomes of host-symbiont interactions, including how Wolbachia spreads to become common.IMPORTANCE Microbes infect a diversity of species, influencing the performance and fitness of their hosts. Maternally transmitted Wolbachia bacteria infect most insects and other arthropods, making these bacteria some of the most common endosymbionts in nature. Despite their global prevalence, it remains mostly unknown how Wolbachia influence host physiology and behavior to proliferate. We demonstrate pervasive effects of Wolbachia on Drosophila temperature preference. Most hosts infected with A-group Wolbachia prefer cooler temperatures, whereas the one host species infected with divergent B-group Wolbachia prefers warmer temperatures, relative to uninfected genotypes. Changes to host temperature preference generally do not alter Wolbachia abundance in host tissues, but for some A-group strains, adult males have increased Wolbachia titer when shifted to a cooler temperature. This suggests that Wolbachia-induced changes to host behavior may promote bacterial replication. Our results help elucidate the impact of endosymbionts on their hosts amid the global Wolbachia pandemic.

Environmental and Genetic Contributions to Imperfect wMel-Like Wolbachia Transmission and Frequency Variation.

Maternally transmitted Wolbachia bacteria infect about half of all insect species. They usually show imperfect maternal transmission and often produce cytoplasmic incompatibility (CI). Irrespective of CI, Wolbachia frequencies tend to increase when rare only if they benefit host fitness. Several Wolbachia, including wMel that infects Drosophila melanogaster, cause weak or no CI and persist at intermediate frequencies. On the island of São Tomé off West Africa, the frequencies of wMel-like Wolbachia infecting Drosophila yakuba (wYak) and Drosophila santomea (wSan) fluctuate, and the contributions of imperfect maternal transmission, fitness effects, and CI to these fluctuations are unknown. We demonstrate spatial variation in wYak frequency and transmission on São Tomé. Concurrent field estimates of imperfect maternal transmission do not predict spatial variation in wYak frequencies, which are highest at high altitudes where maternal transmission is the most imperfect. Genomic and genetic analyses provide little support for D. yakuba effects on wYak transmission. Instead, rearing at cool temperatures reduces wYak titer and increases imperfect transmission to levels observed on São Tomé. Using mathematical models of Wolbachia frequency dynamics and equilibria, we infer that temporally variable imperfect transmission or spatially variable effects on host fitness and reproduction are required to explain wYak frequencies. In contrast, spatially stable wSan frequencies are plausibly explained by imperfect transmission, modest fitness effects, and weak CI. Our results provide insight into causes of wMel-like frequency variation in divergent hosts. Understanding this variation is crucial to explain Wolbachia spread and to improve wMel biocontrol of human disease in transinfected mosquito systems.